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The origin of the helicity hemispheric sign rule reversals in 
the mean-field solar-type dynamo 
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ABSTRACT 

Observations of proxies of the magnetic helicity in the Sun over the past two solar 
cycles revealed reversals of the helicity hemispheric sign rule (negative in the North 
and positive in the South hemispheres). We apply the mean-field solar dynamo model 
to study the reversals of the magnetic helicity sign for the dynamo operating in the 
bulk of the solar convection zone. The evolution of the magnetic helicity is governed 
by the conservation law. We found that the reversal of the sign of the small-scale 
magnetic helicity follows the dynamo wave propagating inside the convection zone. 
Therefore, the spatial patterns of the magnetic helicity reversals reflect the processes 
which contribute to generation and evolution of the large-scale magnetic fields. At 
the surface the patterns of the helicity sign reversals are determined by the magnetic 
helicity boundary conditions at the top of the convection zone. We demonstrate the 
impact of fluctuations in the dynamo parameters and variability in dynamo cycle 
amplitude on the reversals of the magnetic helicity sign rule. The obtained results 
suggest that the magnetic helicity of the large-scale axisymmetric field can be treated 
as an additional observational tracer for the solar dynamo. 

Key words: Turbulence: Mean-field magnetohydrodynamics; Sun: magnetic field; 
Stars: activity - Dynamo 



■ J-i 1 INTRODUCTION 
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Vector magnetographic observations of the solar active re- 
gions show that the distribution of the electric current he- 
licity has a pronounced anti-symmetry with respect to the 
solar equator ([Seehafcr 1990; Pcvtsov et al.|1994 



1995 



fc Zha ng 1998; Kuzany arTet al.||2000| |Hagino fc Sakurai 



Bao 



2005 



Zhang et al. 20101. This phenomenon is called the 



hemispheric sign rule of current helicity. By analysis of the 
photospheric vector magnetograms of active regions it has 
been shown that the current helicity in the northern hemi- 
sphere is mainly negative while in the southern hemisphere 
it is positive. The same hemispheric sign rule was obtained 
from the synoptic mag netic field maps by |Pevtsov fc La-| 
tushko| pOOOt (see, also, |Pevtsov et al.|2001| ). Both kinds of 
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observations deal with the line-of-sight part of current he- 
licity, which can be identified with the total current helicity 
density using the assumption of the spatial isotropy of the 
current helicity distribution. 

It is possible to relate the current helicity density with 
the magnetic helicity density taking into account the the- 
oretical assumption on turbulent nature and isotropy of 
the magnetic fields (see, e.g., Moffatt|1978 Kleeorin & Ro- 
|gachev skii 1999). Magnetic helicity is an integral of motion 
in MHD ( |Woltjer| ( [T958| ; |Moffatt| |l969| ) . This impacts the 
saturation of the magnetic field generation in the large-scale 
helical dynamos (Frisch et al.||1975| Kleeorin & Ruzmaikin 



T982"! IVainshtein fc Kitchatinov||1983| |Kleeorin et al.||200'0 



Brandenburg & Subramanian 2005). Thus, the information 



about the surface distribution of the current helicity density 
and about its evolution with the solar cycle may be impor- 
tant for our understanding of the dynamo processes inside 



the solar convection zone (Kleeorin et al. 2003 Choudhuri 



|et al.||2004 Zhang et al.||2012|. It is also important for un- 



derstanding the processes of the magnetic helicity transport 
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from the convection zone to the corona (Berger & Ruzmaikin 



2000| |Warnecke et al.||2011||Brandenburg et al.|2011 | 

The observations indicate departure from the hemi- 
spheric sign rule ( |Bao et al.|2000||Hagino fc Sakurai|2005[ ). It 
was found that at some periods of the solar cycle the hemi- 
spheric sign rule reverses to the opposite, at least at some 
latitudes and times ( [Zhang et al.|2010| ). It was realized that 
the properties of these reversals may be related with the 
kind of the dynamo operating in the Sun with the distri- 
bution of the dynamo processes inside the convection zone, 
and with the types of the magnetic helicity loss involved in 



the dynamo (see, e.g., Sokoloff et al. 2006 Guerrero et al 



20T0J IMitra et al.||2011[ |Pipin fc Kosov ichcv 20 11b||Zhang 



et al. 20121. These mechanisms do not exclude the local 



processes which may take part in formation of the twisted 
magnetic field at the subsurface layers. Some of them were 
brought attention in the literature and could be considered 
as alternative points of view to the problem (see, e.g., Long 



cope et al 



2007) 



tive points of view to the problem (see, e.g., Long- 
.|1998[ ; Kuzanyan ct al. 2006; Pcvts ovfe Long cope 



The purpose of this paper is to analyze the origin of 
the current helicity sign rule reversals within the framework 
of solar mean-field dynamo models. In our study we exam- 
ine the dynamo distributed over the convection zone. In this 
model the global dynamo wave is shaped by the subsurface 



shear layer (Pipin & Kosovichev 2011c I. Our approach is a 



development of the results of the dynamo model of ( Pipin, 
|Sokoloff, Zhang et al.|[2013[ ) which alleviates catastrophic 
quenching by consideration of total magnetic helicity con- 
servation. We compare our results with that ones for the 
solar dynamo operating in overshoot layer at the bottom of 
the solar convective zone (see Zhang et al.|2012 1 . Our study 
confronts the results of theoretical modeling with available 
observational data from Huairou Solar Observing Station of 
Chinese Academy of Sciences. 



2 BASIC EQUATIONS 

The details of the model can be found in our previous pa 
pers (see, e.g., |Pipin et al.|2012|[Pipin|2013| and |Pipin|2008 



hereafter P08). Here, we briefly outline the basic framework. 
We study the mean-field induction equation: 



dB 

dt 



= V x (£ + UxB) 



(1) 



where U is the mean velocity (differential rotation); B is the 
axisymmetric magnetic field: 



B = e^B + V x 



r sin 8 ' 



(2) 



where 8 is a polar angle and r is a radial distance; £ — u X b 
is the mean electromotive force, with u and b are being the 
fluctuating velocity and magnetic field, respectively. Using 



the mean-field magnetohydrodynamic framework ( Krause & 
Radler 1980) we write the £ as follows: 



(a lJ+1 ^\B 



J Bj - \r)i ih + rjijl) 



//,< +ri$'k)*,B... (3) 

where the turbulent kinetic coefficients are: the a effect, am; 
the turbulent pumping j>. '; the anisotropic diffusivity rjijk 
and the S dynamo effect (Radler 19691, r)^A. They depend 



on the parameters of the turbulent convection, like the mean 
density and turbulent diffusivity stratification, the Corio- 
lis number fl* — 2r c flo, where r c is the typical convective 
turnover time, and f2o is the global angular velocity. 

The a effect includes the hydrodynamic and magnetic 
helicity contributions, 

(4) 



a 



n ■ 2 Q (H) , (Al) 

v, = C a sin 8a.). + a) 



»J 



»J 



where the hydrodynamic part of the a-effect is defined by 

(H) ' 

\ - ■ 
13 

a). , 7}. , rjijk and rf^A are given in Appendix. The con- 



»j 



tribution of small-scale magnetic helicity \ — a • b (a is the 
fluctuating magnetic vector-potential) to the a-effect is de- 
fined as follows (see, P08): 



AM) 



= 2/ a w * 



(a); 



XTc 



2f[ a) e,e. 



XTc 



(5) 



"HOPP J1 Vo^ 2 

The principal nonlinear feedback of the large-scale magnetic 
field to the a-effect is due to a dynamical quenching because 
of the generation of the magnetic helicity by the dynamo 
(jFrisch et al.||1975| |Kleeorin fc Ruzmaikin||1982| [Branden- 
burg fc Subramanian|2005 1 . The relation of magnetic helic- 
ity on the large- and small scales ( Hubbard & Brandenburg 
2012; |Pipin|2013| is governed by the equation 

9x d(A-B) 

dt 



i)t 



R„ 



r/B J V-T x 



(6) 



where JF X = — r\ x V (x + A ■ B) is the diffusive flux of the 



total magnetic helicity ( [Mitra et al. 20101, and rj x is the 
turbulent diffusion coefficient for the magnetic helicity. In 
the paper we use R m = 10 6 and i] x = O.lr^T, where rjr is 
the turbulent diffusivity profile (see, Appendix). For the ax- 
isymmetric magnetic fields the large-scale magnetic vector- 
potential is 

e (? , 



A = e T + yP 



-A + re r P. 



(7) 



rsin.8 

The toroidal part of the vector potential is governed by the 
dynamo equations. The poloidal part of the vector potential 
can be restored from equation V x (rP) = e^B. We matched 
the potential field outside and the perfect conductivity at the 
bottom boundary with the standard boundary conditions. 
For the magnetic helicity we employ \ — at the bottom of 
the convection zone. The paper elaborates two kind of the 
surface boundary conditions for the magnetic helicity: 

r; x V r (x + A-B)| r=r = 0, (8) 



rj x \7 r 



(pP) 



= 0. 



(9) 



We call the model that satisfies the boundary conditions 
Eq. pll as the model Bl, and similar, the model B2 is re- 
ferred to the Eq. pi. We set the seed magnetic field of 
the preferred dipole parity and with small admixture of the 
quadrupole one to check the parity preference when the so- 
lution reaches the steady state. 

The construction of the radial profiles for the turbulent 
coefficients, which are involved in the mean electromotive 
force, remains rather arbitrary for various kinds of the dy- 
namo models. In our models we use the solar convection 
zone model computed by |Stix| ( |2002| ) . In the paper we use 
the same profiles for the turbulent coefficients as in our pre- 
vious papers (see, Pipin|2013 and Appendix therein). 
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Figure 1. Snapshots of the magnetic field and magnetic helicity evolution inside the convection zone: Top panel shows the field lines of 
the poloidal component of the mean magnetic field and the toroidal magnetic field (varies ±lkG) is shown by color; Bottom panel shows 
the small-scale magnetic helicity density (contours) and the large-scale magnetic helicity density (color). Both quantities vary with the 
same magnitude. 








Figure 2. The left column shows the results for the model Bl. Panel (a) shows the time-latitude diagram for the current helicity 
(background image). Panel (b) shows the toroidal magnetic field variations at r = 0.95-R. panel (c) shows variations of the small-scale 
magnetic helicity and the toroidal magnetic field inside the convection zone at the latitude 30°. Panels (d, e, f) show the same results 
for the model B2. 
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3 RESULTS 

Fig. IT] shows the snapshots of the magnetic field and mag- 
netic helicity evolution in the North hemisphere for the 
model Bl which uses Eq. dSu - The qualitatively similar re- 
sults can be obtained for the model B2. Here we see, that 
the spatial patterns of the small-scale magnetic field follow 
the evolution of the large-scale magnetic helicity and the lat- 
ter propagates with the toroidal part of the dynamo wave 
from the bottom of the convection zone to the surface. The 
dynamo wave has the equatorial and the polar branches. 
Near the surface the equatorial branch dominates. The hemi- 
spheric helicity rule suggest that the small-scale helicity is 
negative at the North and positive at the South hemisphere. 
Fig. [T] shows that in the upper part of the convection zone 
the helicity rule is valid in the most phases of the cycle. 
In the upper part of the convection zone the reversal sign 
of the small-scale magnetic helicity regions appears at the 
high latitudes when the dynamo wave of the toroidal mag- 
netic field comes to the subsurface shear layer. At the equa- 
torial latitudes the reversal sign of \ occurs at the decaying 
phases of the dynamo wave cycle. One can see that the signs 
of the large and small-scale helicities are spatially related. 
It is shown in Fig. [2] that shows the time-latitude and the 
time-radius variations of the magnetic field and magnetic 
helicity near the surface. The Figure also demonstrates the 
effect of the boundary condition change for the magnetic he- 
licity. For the boundary condition Eq. (TSj) the regions with 
reversed sign of the small-scale magnetic helicity penetrate 
into the surface while the condition Eq. |9| quenches this 
penetration. We find that the patterns of the reversed sign 
of the magnetic helicity are located at the edges of the but- 
terfly wings of the time-latitude diagrams for the large-scale 
toroidal magnetic field. The novel feature which is demon- 
strated by the Fig. [2] is the time-latitude diagram for the 
large-scale magnetic helicity which is attributed to the ax- 
isymmetric magnetic field. It is seen that within the current 
model its distribution is closely connected with the distri- 
bution of the small-scale helicity. It is believed that the cur- 
rent helicity of the surface magnetic field is the observational 
proxy for the magnetic helicity \. We note that our model 
uses the full information about the large-scale magnetic he- 
licity. 

The difference in penetration of the magnetic helicity 
to the surface results in difference in the distribution of the 
effective a effect near the surface. This issue is recently dis- 



cussed by Kapyla, et al 



Pipin, Sokoloff, Zhang et 



(|2012| ) 

al. (20131 and Pipin (20131. Fig. 13] shows the snapshots of 
the Q, i<? , and the small-scale magnetic helicity profiles for the 
different phase of the cycle at the latitude 45°. We find that 
for the model Bl the a-effect can be negative at the cer- 
tain phases of the cycle and it has the sharp positive profile 
near the surface. The model B2 has the negative a-effect for 
r > 0.92R with the sharp negative profile near the surface. 
The abrupt growth of the a effect amplitude near the sur- 
face is because of the factor (p£j in the definition, see 
Eq. |5| . It remains the matter of the direct numerical simula- 
tions to justify the correct choice of the boundary condition 
for the magnetic helicity. The model B2 has the zero bound- 
ary condition for the derivative of the small-scale current 
helicity at the top, see Eq. [9] It is found that for the con- 
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Figure 3. Panels (a,b) show variations of the a-effect and the 
small-scale magnetic helicity at the latitude 45° for the model Bl 
(the Eq. (ISp) , and the panels (c,d) show the same for the model 
B2 (the Eq. {§}). 



dition V r xlr — r ~ ® ^he negative part of the alpha-effect 
near the surface is stronger than one in model B2. 



3.1 Impact of dynamo fluctuations on the helicity 
patterns 

The sign reversals of the helicity rule can be due to random 
fluctuations in the dynamo parameters and due to some ran- 
dom processes which generate the magnetic helicity indepen- 
dent of the large-scale dynamo. In this subsection we exam- 
ine the effect of fluctuations in the dynamo parameters on 
the magnetic helicity distribution variations. We exploit here 
a scenario (Moss et al.|2008 Usoskin et al.|2009 Pipin et al 



|2012[ ) with fluctuations of a-effect as a possible source of the 
solar activity cycle parameters from one cycle to another. We 
introduce random non-symmetric about equator variations 
of the a-effect, C a = C a (1 + 0.2 (£y6 (n) + £<?6 (-ft))), 
where \x = cos#, is the Heaviside function and |£s,jv| < 
2o"(£s,at) is the random Gaussian noise with the randomly 
floating phase and with the mean memory time equals to 
the dynamo cycle length. In this subsection the model Bl 
is discussed as it shows the stronger reversals of the helicity 
rule than the model B2. 

We found that the reversals of the helicity rule are 
stronger during the periods of the grand minimum which 
are also related to the periods of the strong hemispheric 
asymmetry in the magnetic activity. Fig. HI shows variations 
of the integral parameters of the model for the near sur- 
face magnetic field. In our results we show the parity index, 
that determines the symmetry of the toroidal magnetic field 
about equator, with the value —1 corresponds to the dipo- 
lar symmetry of the near surface toroidal magnetic fields 
and the value 1 corresponds to the quadrupolar symmetry. 
The sunspot number was simulated in following to |Pipin 



et al. ( 2012 1. We also show the integral magnetic helicity for 



each hemisphere. The magnitude of the helicity variations 
is in agreement with the observational constraints obtained 
by |Berger fc Ruzma ikin (2000). Variations of the magnetic 
helicity go in anti-phase at the large and small scales, be- 
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Figure 4. The integral characteristics of the magnetic activity near the surface, r = 0.95i? : a) the parity index determining the symmetry 
of the toroidal magnetic field about equator, -1 corresponds to the dipolar symmetry and 1 to the quadrupolar; b) the simulated sunspot 
number (SN); c) the latitudinal integral of the large-scale magnetic helicity, dashed line the North hemisphere, the dash-dotted line the 
South one; d) the same as (c) for the small-scale magnetic helicity; e) the same as (b), (the SN), the solid line is the total SN, the dashed 
line is the SN for the North hemisphere and the dash-dotted line is the SN for the South hemispere. 



cause it is prescribed by Eq.™. Nevertheless, for each hemi- 
sphere, there is a difference between the evolution of AB 
and x = a 'b- The small-scale helicity, x does change the 
sign in a course of the solar cycle and the large-scale helic- 
ity A-B almost does not. This is similar to results shown in 
Fig. [2] where we see that reversals of the helicity rule is much 
stronger for the small-scale helicity than for the large-scale 
one. Another interesting results is that the maxima of the 



integral large-scale magnetic helicity are approximately cor- 
responded to the maxima of the decay rate in the simulated 
sunspot activity (cf, Figs. En, d and Fig. He). This is due of 
the oscillatory character of the dynamo and delay between 
the activity of the major components of the large-scale mag- 
netic fields which are related to the toroidal magnetic field 
and the large-scale toroidal vector-potential determining the 
poloidal magnetic field. 
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Finally, Fig. [5] shows comparison of the results for the 
simulated time-latitude diagrams for the toroidal magnetic 
field and the magnetic helicity with results of the current 
helicity observations reported by Zhang et al. (20101. We 



used systematic series of vector magnetographic observa- 
tions of solar active regions by 35 cm filter type SMTF 
telescope at Huairou Solar Observing Station of Chinese 
Academy of Sciences. The data set comprises 6205 individ- 
ual magnetograms of active regions more or less homoge- 
neously covering the 18 year period of 1988-2005, which is 
almost two sunspot cycles. The data have been grouped and 
averaged into statistically significant sub-samples in time- 
latitudinal bins (2 years in time and 7 degrees in helio- 
latitude), see Zhang et al. (20101 for details. We have sub- 



sequently smoothed the data by using standard IDL linear 
interlopation for retaining only global features of the time- 
latitudinal distribution of helicity. 

The results of the dynamo model are shown for the pe- 
riod of the grand minimum. It is the same period as dis- 
cussed for the Fig. HTc,d,e) above. The simulated butterfly 
diagrams are in visible qualitative agreement with the ob- 
servations. One can see that the model keeps the basic anti- 
symmetry of helicity (negative in the North and positive 
in the South, i.e. the so-called hemispheric sign rule), how- 
ever with evolution it shows various deviations from perfect 
periodicity (e.g., longer cycles, suppression of activity, asym- 
metry in the phases of growth and decay, asymmetry in the 
shape of wings on butterfly diagrams etc). It looks plausible 
that long and weak cycles are associated with larger areas 
of helicity of the sign opposite to the hemispheric sign rule. 



4 DISCUSSION AND CONCLUSIONS 

The available bulk of the current helicity data covers two ac- 
tivity cycles and the transition to the following activity cycle 
which has been quite unusual. The observed helicity butter- 
fly diagrams demonstrate that the size of the areas with the 
opposite helicity signs in the later cycle differs substantially 
from the first one. We studied the origins of the reversals of 
the magnetic helicity sign in the mean-field solar dynamo. 
The evolution of the magnetic helicity in the model subject 
to the global constraint of the magnetic helicity conserva- 
tion law. The nonlinear feedback of the large-scale magnetic 
field to the a-effect is described by dynamical quenching 
due to the constraint of magnetic helicity conservation. The 
magnetic helicity, %i ls subjected to the conservation law. 

In the model, the sign reversals of the small-scale mag- 
netic helicity are always related with the sign reversals of 
the large-scale magnetic helicity. This is due to the mag- 
netic helicity conservation constraint. The result develops 



the simple model by Xu et al. (2009). The idea was recently 



elaborated by Zhang et al. (20121 for the toroidal part of 



the current helicity. Our model employs the total large-scale 
magnetic helicity and not only its toroidal part. Taking into 
account Eqs.(2|7l we get the large-scale magnetic helicity 
formula for the spherical coordinates: 



A B 



where, B = B 



(10) 



AB P dA 

rsin8 rsin# 88 ' 

4, = - — - and B r = . -— . We have to 
08 r 1 sin 8 08 



notice that the magnetic helicity of the large-scale axisym- 
metric field can be restored from observational tracers of 
B4 and B 



Seehafer 



either from the vector magnetograms 
(1990) or from the line-of-sight magnetic observations, e.g., 
using method by |Pevtsov et al. (2001 1. The toroidal part of 
the potential can be restored from the surface distribution of 
the B r as A (8) = J r 2 sm8B r dd which is equivalent to the 
flux going outside of the Sun, and, similar, we can restore 
poloidal part of vector-potential using P — — J B^dd. Note, 
that the total helisity remains zero because of the equatorial 
symmetry of the axisymetric magnetic field. Hence, the us- 
ing of the Coulomb gauge is justified and the determination 
of the latitudinal distribution of the magnetic helicity of the 
axisymmetric large-scale magnetic fiel is unique. 

Therefore, the observations can give an information 
about the magnetic helicity of the large-scale magnetic fields 
of the Sun. Our results indicate (see, Figs. B[c,d,e) that the 
reversal of magnetic helicity and lower values of integral he- 
licity may proceed the lower amplitude of cyclic dynamo 
activity. 

We found that in the models the wave of the reversed 
magnetic helicity sign propagates from the bottom of the 
convection zone. Similar property was recently found by 



Warnecke et al. (20111 in direct numerical simulations. The 



models Bl and B2 show the possibility as for the strongly 
positive as well for the negative dynamical a-effect near the 
top of the convection zone. The model employs the sub- 
surface rotational shear having the negative radial gradient 
of the angular velocity. Therefore, in case Bl the dynamo 
wave penetrates closer to equator than the model B2 be- 



cause of the Parker- Yoshimura rule ( Parker 1955 Yoshimura 



1975[). The numerical sim ulations demonstrate a similar ef- 
fect ( |Kapylaet al.|2012 |. 

Generally, the sign reversals of the magnetic helicity 
are stronger in the model Bl than in the model B2. These 
two models have different boundary conditions for the mag- 
netic helicity at the top. In the model Bl the diffusive flux 
of the large-scale helicity from the surface is balanced by 
a counterpart from the small-scale helicity. Therefore, inte- 
grating the Eq. (18} from some level r = ro the top r = r e 
we have x e = Xo + -^-o ■ Bo, where we use the magnetic 
boundary conditions as well. Thus, in the model Bl, the 
boundary conditions support the penetration of the local 
helicity \o + Ao • Bo (governed by Eq. p| from depth to 
the surface. For the boundary condition Eq. (T9J) we have 
V r ftc = at the top. This is the same as V r X = X-Mp^ 2 ), 
where A(p£ 2 ) = V r log(p^ 2 ). Thus, at the near surface level, 
in the model B2 the small-scale magnetic helicity is deter- 
mined by the profile of A(p£ 2 ) and not by Eq.™. The further 
study requires clarification of the issue if the boundary con- 
ditions impact the sign reversals of the magnetic helicity. 
We can make conjecture that the change in the boundary 
conditions results in the larger or smaller sign reversals of 
the magnetic helicity at the surface. 

The main results of the paper can be summarized as fol- 
lows. The current model suggests that the reversal of the sign 
of the small-scale magnetic helicity follows the dynamo wave 
propagating inside the convection zone. This was also sug- 
gested by the numerical simulations ( Warnecke et al.|2011 1 . 
Therefore, the spatial patterns of the magnetic helicity re- 
versals reflect the processes which contribute to generation 
and evolution of the large-scale magnetic fields. At the sur- 
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Figure 5. The time-latitude diagram for the current helicity as inferred from the solar active region observations at Huairou Solar 
Observing Station. Sunspot density which traces toroidal magnetic field is shown by colors while contour lines show current helicity, the 
vertical color bar on the right side scales the magnitude. 



face the patterns of the helicity rule reversals are determined 
by the magnetic helicity boundary conditions at the top of 
the convection zone. The model suggests that the magnetic 
helicity of the large-scale axisymmetric field can be used as 
an additional observational tracer for the solar dynamo. 
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APPENDIX 

Here we describe some details of the dynamo model that can 
be al so fond in ( |Pipin et a,l7||2012| |Pipin||2013| ) and ( |Pipin 
2008 1 (hereafter P08). The hydrodynamic part of the tensor 
Qij is represented by a\, (P08): 



(H) 
M ■ ■ 
<-3 



+ 
+ 



^{3^(/W(e.AW)+/W(e.AW))}( 

{ 3?7T (/f)(-A W )+/i a) (-A (u) ))} 
3^(eiAW+c-AM)/e (B) 
(«A« + e,AM)/M}, 



where, A^ p ' = Vlogp is the inverse density stratification 
height, A> u ' = |Vlog [rj^. ) is the same for the turbulent 
diffusivity, e = n/|J7| is a unit vector along the axis of 
rotation. The turbulent pumping, 7} , depends on mean 
density and turbulent diffusivity stratification, and on the 
Coriolis number Q,* = 2r c Qo where r c is the typical con- 
vective turnover time and ilo is the global angular velocity. 
Following the results of P08, 7}- is expressed as follows: 

7^ = 3 VT {ft ) Ai P) + f[ a) (e-A^)e n }s lnj (12) 

- Wf^'pf. p A (p) 

- 3vt (e - 1) {/ 2 W A^ + /<"> (e ■ A'"') e„} e inj . 

The effect of turbulent diffusivity, which is anisotropic 
due to the Coriolis force, is given by: 

Vijk - 3vt I f 2/ 1 (a) - / 2 (d) J e ijk - 2f[ a ' aenSnjk > ■ (13) 

We also include the nonlinear generation effects which is in- 
duced by the large-scale current and the global rotation that 
is usually called as the Q x J effect or the 5 dynamo effect 
( Radl er[l969 1. It is supported by the numerical simulations 
JKapyla et al.|2008l |Schrinner||2011[ ). P08 suggested that: 



(S) n ,-, r(d) J ~(™)r , ~(w)Bi.Bk 

%k = ZnrCsU 'ej ^ <p\ '5 ik + <p\ '—s r 



(14) 



where, Cs measures the amplitude of the Q x J effect, 
(p 2 7 (/3) are normalized versions of the magnetic quench- 



,(») 



ing functions tp 2 7 given in P08. They are defined as follows, 



„(») 



£27 W) = 1^2 7 (P)- The last term in Eq.(14l is the non- 



linear contribution to the fl x J-effect. Its structure is the 
same as for the a effect because the associated electromo- 
tive force is proportional to I^tC^/I (p™ Bi ( e ' V) logB 
(see details in P08). Thus this effect works similar to the a 
effect that is excited by the nonlinear buoyant instabilit y of 
large-scale magnetic field. The functions /Sj^^iin Eqs( 13D 
depend on the Coriolis number. They can be found in PD8 
(see also, |Pipin fc Kosovichev|2011b| |Pipin fc Sokoloff|2011 1 . 

The mixing-length is defined as £ = «mlt A™ , 
where A^ p ' — V logp is the inverse pressure variation height, 
and q:mlt = 2. The turbulent diffusivity is parameterized in 



the form, r\T 



^vVt ' 



where rfj, 



is the charac- 



3/ D „ (r) 

teristic mixing-length turbulent diffusivity, £ is the typical 
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Helicity hemispheric sign rule 9 



correlation length of the turbulence, C n is a constant to con- 
trol the efficiency of large-scale magnetic field dragging by 
the turbulent flow. Also, we modify the mixing-length tur- 
bulent diffusivity by factor f ov (r) = f + exp (50 (r ov — r)), 
r ov — O.725-R0 to get the saturation of the turbulent pa- 
rameters to the bottom of the convection zone. The latter is 



suggested by the numerical simulations (see, e.g., Ossendri 



jver et al. 2001 2002 Kapyla et al. 20081. The results do 



not change very much if we apply A^ u ' = CV,Vlog I rf^. ) 

with C v ^ 0.5. For the greater C v we get the steady non- 
oscillating dynamo which is concentrated to the bottom of 
the convection zone. The purpose to introduce the additional 
parameters like C v = 0.5 and fov(r) is to get the distribution 
of the a effect closer to the result obtained in the numerical 
simulations. 

The bottom of the integration domain is rt — 0.715-Rq 
and the top of the integration domain is r e = O.99i?0. 
The choice of parameters in the dynamo is justified by our 
previous studies ( Pipin fc Kosovichev||2011a l, where it has 
been shown that solar-types dynamos can be obtained for 
C a /C$ ^ 2. In those papers we find an approximate thresh- 
old C a ~ 0.03 for a given value of diffusivity dilution factor 
C v — 0.05. The latter was chosen to tune the solar cycle 
period. 
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